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Using density functional theory, we determine parameters of tight-binding Hamiltonians for a 
variety of Fabre charge transfer salts, focusing in particular on the effects of temperature and 
pressure. Besides relying on previously published crystal structures, we experimentally determine 
two new sets of structures; (TMTTF) 2 SbF 6 at different temperatures and (TMTTF) 2 PF 6 at various 
pressures. We find that a few trends in the electronic behavior can be connected to the complex 
phase diagram shown by these materials. Decreasing temperature and increasing pressure cause the 
systems to become more two-dimensional. We analyze the importance of correlations by considering 
an extended Hubbard model parameterized using Wannier orbital overlaps and show that while 
charge order is strongly activated by the inter-site Coulomb interaction, the magnetic order is only 
weakly enhanced. Both orders are suppressed when the effective pressure is increased. 



I. INTRODUCTION 

Quasi-one dimensional organic salts formed from 
tetramethyltetrathiafulvalene (TMTTF) molecules - also 
known as Fabre charge-transfer (CT) salts - have been in- 
tensively investigated in the last two decades since they 
exhibit a rich variety of phases like antiferromagnetism, 
superconductivity, charge ordering, spin-density wave or- 
dering or spin-Peierls behaviorpJ' 3 -' Such phases can be 
driven both by external (physical) pressure as well as 
by chemical pressure (see Fig. [T]). Even though a few 
successful models have been proposed for the description 
of these systems, constructing a consistent microscopic 
picture of the relationships between the various phases 
remains a challenge 

The primary avenue of the present work is to under- 
stand the microscopic origin of the close competition be- 
tween the different phases in these compounds as a func- 
tion of chemical and external pressure as well as tempera- 
ture. For tha t, w e performed ab initio density functional 
theory (DFTpED anc j model Hamiltonian calculations for 
several Fabre CT salts whose crystal structures were de- 
termined at different temperatures and pressure and in- 
vestigated variations of their electronic properties with 
temperature and pressure!^ By computing the real-space 
overlaps of Wannier orbitals for the bands near the Fermi 
level, we parametrize a two-band tight-binding Hamilto- 
nian model for the various systems and examine the dif- 
ferences in their electronic hopping parameters. In this 
way, we can connect structural and chemical modifica- 
tions with changes in the electronic properties. Further- 
more, in order to analyze some of the preferred orderings 
we consider a description of the Fabre CT salts in terms of 
an extended Hubbard Hamiltonian including on-site and 
inter-site Coulomb interaction terms. The kinetic part 
of this model is given by the computed hopping parame- 
ters for the various compounds. We discuss spin-spin and 
charge-charge correlation properties by diagonalizing the 



model. 

The work is organized as follows. Sections [H] and |III| 
are dedicated to the description of the computational de- 
tails as well as the crystal structure of the Fabre CT salts. 
In Sections IV V]and VI we present, respectively, our re- 
sults on the electronic structure, tight-binding models as 
well as extended Hubbard models for a few members of 
Fabre CT salts. Discussion and conclusions are given in 
Sections EED and EHO 



II. COMPUTATIONAL DETAILS 

The electronic structure calculations presented here 
were performed in an all-electron full-potential local or- 
bital basis using the FPLO package.^ The densities were 
converged on an (8 x 8 x 8) k mesh, using a generalized 
gradient approximation (GGA) functional.^ 

For materials without published hydrogen coordinates, 
hydrogen atoms were placed according to the expected 
bond lengths and angles of a methyl group (C-H distance 
1.1 A, C-C-H angle 109°). With the bond length and 
angle fixed, one has the freedom to choose the rotation 
angle of the set of hydrogen atoms on each methyl group 
around the C-C bond. We chose this angle such that 
one hydrogen is as far out of the plane of the molecule 
as possible. We tested the effect of this choice on the 
band structure and found no contribution to the bands 
of interest near the Fermi level since in this energy region 
only TMTTF bands of n origin are involved. 

Some of the structures had suspicious bond lengths 
and angles, therefore the atomic coordinates were relaxed 
with DFT using the Vi enna, ab initio Simulation Package 
(VASP, ve rsion 5.2.11) EM3 with a projector-augmented 
wave basis liSliU We used the GGA functional,^ and in- 
cluded Van der Waals corrections 2 ^ for the relaxations. 
We performed two kinds of relaxations; in one relaxation 
we kept sulfur and the heavy anion atoms coordinates 
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FIG. 1. Temperature-pressure phase diagram for the TMTTF 
and TMTSF charge transfer salts, as first suggested in Ref. 
[H and refined by many others. Position in the phase dia- 
gram can be tuned by physical pressure, or chemical pressure 
(changing anion) . The ambient pressure position for each salt 
is indicated with an arrow above the diagram. An increase of 
pressure (external or chemical), causes the system to be less 
one dimensional. The position of (TMTTF) 2 C10 4 is not well 
known. In the phase diagram its approximate position has 
been indicated with a dashed arrow. The possible phases are 
charge ordered (CO), Mott insulating (MI), antiferromagnetic 
(AF), spin Peierls (SP), spin density wave (SDW), supercon- 
ducting (SC) and ID, 2D or 3D metal. 




FIG. 2. Crystal structure of the Fabre CT salts projected into 
the ac plane. The organic molecules form 7r-stacked ID chains 
along the crystal a direction (with a slight zig-zag pattern), 
and form layers parallel to the ab plane. These organic layers 
alternate with anion layers (with the anions centered on the 
pink As sites) stacked in the c direction. Grey atoms are 
carbon, yellow are sulfur, while hydrogen atoms are shown in 
white. 



III. CRYSTAL STRUCTURE 



fixed, and in the other relaxation all atom positions were 
relaxed. The differences between the two relaxed struc- 
tures were minimal. The atomic coordinates were con- 
verged to an energy difference of 1 meV on a (5 x 5 x 5) 
/c-mesh, with a plane wave cutoff energy of 500 eV. 

The tight-binding parameters were obtained by con- 
structing Wannier orbitals for the TMTTF bands at the 
Fermi level and computing real-space overlaps, as imple- 
mented in FPLO. Another way to generate these parame- 
ters is to fit the band structure of the model Hamiltonian 
to the DFT bands. The latter method can become dif- 
ficult when many hopping parameters need to be fitted; 
there can be a number of solutions which reproduce the 
DFT bands equally well, but differ in physical details 
(such as relative strengths of certain bonds). By using 
Wannier orbital overlaps we can be sure that our param- 
eter values have a clear physical interpretation. 

The exact diagonalization of the extended Hubbard 
model was performed by considering system sizes of 4 x 4 
TMTTF sites with periodic boundary conditions. 



The Fabre CT salts consist of alternate layers of 
TMTTF molecules (cations) and monovalent anions, 
stacked in the c direction (see Fig. |2j). In between the 
cation layers, the planar TMTTF molecules form ir- 
stacked one dimensional chains in the a direction with 
a slight 'zig-zag' arrangement. There is a charge trans- 
fer of one electron from each (TMTTF) 2 dimer to each 
anion, i.e. the TMTTF molecules carry half a hole on 
average. There are two classes of anions: those that con- 
form to the P 1 symmetry of the TMTTF part of the 
crystal (such as PFg), and those that break that sym- 
metry (such as CIO4). Also, the anion species influence 
the proportions of the unit cell, as well as the intra- and 
inter-chain hopping strengths. The inter-chain hopping 
strengths are not only determined by the distance be- 
tween the TMTTF molecules, but also by changes in their 
zig-zag arrangement (that is to say, how far away from a 
perfectly aligned stack they are, and in what direction). 

TMTTF molecules within a chain show a slight dimer- 
ization along the chain. We can quantify this structural 
dimerization as the difference between the larger dimer- 
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ization distance of adjacent TMTTF molecules, d±, and 
the shorter dimerization distance, doi normalized by the 
sum of the two distances: 

a _ di — dp 

O struc i , 7 i \ ) 

di + d 

These distances are defined as the distances between the 
centres of mass of the C and S atoms in each TMTTF 
molecule. 

Table|T]shows the structural dimerization of the materi- 
als investigated in this work. This table also includes the 
electronic dimerization, which is introduced in Sec. V. 
In general, the structural dimerization increases slightly 
with increasing temperature. 

Cooling from room temperature to T = 4 K, 
(TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 both show charge 
ordering phase transitions and spin-Peierls transitions. 
These ordering transitions are not visible in the crys- 
tal structures; there are no significant changes in the 
structural dimerization from T= 300 K to T= 4 K. 
(TMTTF) 2 SbFg shows similar changes in the structural 
dimerization between 100 K and room temperature, and 
does not go through any ordering transitions in this 
range. 

IV. ELECTRONIC STRUCTURE 

In the following we examine the electronic properties 
of (TMTTF) 2 PF 6 in detail and will use this analysis as a 
baseline for understanding the Fabre CT salts. In Fig. [3] 
(a) we present the band structure and density of states 
of (TMTTF) 2 PF 6 in a window of energy [-4eV,4eV] 
around the Fermi level. The bands have been drawn 
along the high symmetry path shown in the Brillouin 
zone in Fig. [3] (b). It is clear from the partial density 
of states that near the Fermi level all the bands are pre- 
dominantly due to the TMTTF molecules. In fact, the 
nearest anion bands are about 4.1 eV below the Fermi 
level, and more than 10 eV above it. The two 3/4-fillcd 
organic bands near the Fermi energy are a common fea- 
ture of the Fabre salts as a result of hole doped pairs of 
TMTTF molecules. In (TMTTF) 2 PF 6 these bands are 
well separated from the rest of the bands, with gaps of 
more than 1 eV to the lower valence bands and upper 
conduction bands respectively. The size of the gaps vary 
with anion type, and sometimes anion bands cross the 
two TMTTF bands (as in the case of (TMTTF) 2 Br). 

The quasi-one dimensionality of this system is mani- 
fested in the band structure (Fig. [3] (a)) where we find 
very little dispersion in the k y and k z directions and 
bands only cross the Fermi level in the k x direction. This 
can be also observed in the Fermi surface cut at k z = 
shown in Fig.[3](c). This quasi-one dimensional behavior 
is a typical feature of the Fabre CT salts. 

In order to further characterize the electronic struc- 
ture of these systems, we generate Wannier orbitals for 
the two organic bands near the Fermi level as described 
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Ref. 


SbF 6 (100 K) 


0.007 


0.042 


new 


SbF 6 (140 K, sample 1) 


0.011 


0.067 


new 


SbF 6 (140 K, sample 2) 


0.013 


0.094 


new 


SbF 6 (180 K) 


0.020 


0.115 


new 


SbF 6 (200 K) 


0.023 


0.141 


new 


SbF 6 (300 K, sample 1) 


0.047 


0.279 


new 


SbF 6 (300 K, sample 2) 


0.041 


0.298 


new 


AsF 6 (4K) 


0.007 


0.100 


no 


PF 6 (4K) 


0.009 


0.126 


122 


AsF 6 (300 K) 


0.041 


0.110 I22T 


PF 6 (300 K) 


0.040 


0.230 


new 


PF 6 (300 K, 0.3 GPa) 


0.018 


0.577 


new 


PF 6 (300 K, 0.6 GPa) 


0.016 


0.595 


new 


PF 6 (300 K, 0.9 GPa) 


0.002 


0.477 


new 


PF 6 (300 K, 1.5 GPa) 


0.003 


-0.454 


new 


PF 6 (300 K, 2.0 GPa) 


0.010 


-0.397 


new 


PF 6 (300 K, 2.7 GPa) 


0.024 


-0.183 


new 


Br (300 K) 


0.019 


-0.189 HSF 


C10 4 (300 K) 


0.040 


o.6i6 m* 


BF 4 (100 K) 


0.020 


-0.054 


ma 


BF 4 (300 K) 


0.028 


0.336 


EH 



TABLE I. Structural and electronic dimerization of the Fabre 
CT salts considered in the present work. The structural 
dimerization of the TMTTF molecules is defined in Eq. |TJ, 
and the electronic dimerization is defined in Eq. Q in Sec. 
V. References marked with * have no (or unrealistic) pub- 
lished hydrogen coordinates. Note that CIO4 and BF4 are 
tetrahedral anions, and so do not conform to the reported 
Pi symmetry (the anions do not have the required inver- 
sion symmetry). The anion ordering transition only occurs in 
these systems where the anion does not have inversion sym- 
metry. The change in sign of <9 e i ec indicates that the shorter 
bond has the smaller t. 

in Section [l| An example is shown in Fig. [4j These 
bands have the symmetry of the TMTTF highest oc- 
cupied molecular orbitals (HOMOs), partially depopu- 
lated by the charge transfer of one electron from a pair of 
TMTTF molecules to the anion layer. These two bands 
determine the low-energy physics of these systems and, 
in what follows, we shall concentrate on the analysis of 
this band manifold. We note that we are not considering 
DFT calculations beyond GGA and therefore leave cor- 
relation effects (beyond GGA) to be explicitely treated 
in the model calculations. 



V. TIGHT-BINDING MODEL 

Wannier orbitals form a natural basis for a tight- 
binding model. By computing overlaps between the or- 
bitals, we can parameterize the two HOMO bands at the 
Fermi energy in terms of a two-site tight-binding Hamil- 
tonian where the lattice sites arc defined as the centers 
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of mass of the two TMTTF molecules in each unit cell: 



H 



N 



(•»,i>w 



(2) 



/i is the on-site energy, are hopping parameters be- 
tween sites i and j and the sum over indicates 
that only hoppings up to the N th nearest neighbor are 
included. In listing hopping parameters, we will use 



U 



t a { r ij) 



(3) 



where distances between TMTTF centers of mass 

and a — 0,1,2,... counts neighbour distances in as- 
cending order. In the discussion that follows, we in- 
clude hoppings up to the 8 th nearest neighbor (N = 8). 
These 8 hopping terms do not include any inter-layer 
hopping, and therefore the resulting tight-binding bands 
have no dispersion in the k z direction. The resulting 
tight-binding parameters for the eight shortest inter-site 
distances are shown in Table [TXJ The longer hopping 
terms are of the same order as ti or smaller. 



FIG. 3. Electronic properties of (TMTTF) 2 PF 6 (T = 4 K 
structure), a) Band structure and density of states, b) Path 
through fc-space considered for the band structure plotting, 
c) Fermi surface in the k z = plane . The total density 
of states is shown in black and the partial density of states 
of the anions (increased by a factor of 100) is shown in red 
(dashed) . The partial density of states shows that within this 
energy window, all of the bands in this energy window have 
predominantly TMTTF character, and the two bands at the 
Fermi level are nearly purely TMTTF. 




FIG. 4. Wannier orbital for the (TMTTF) 2 PF 6 bands near 
the Fermi level. It is clear that this Wannier orbital has the 
symmetry of the HOMO of a TMTTF molecule in the gas 
phase. The other Wannier orbital needed to describe the 
two organic bands corresponds to the HOMO of the second 
TMTTF molecule in the unit cell (and is related to the first 
one by inversion symmetry). This is in agreement with the 
information in the partial density of states (Figure [3] (a)); the 
orbitals near the Fermi energy are predominantly of TMTTF 
nature. 



A. Anion dependence of the structural and 
electronic properties 

In Figure [5] we show the band structure of the vari- 
ous Fabre CT salts considered in this study with crys- 
tal structures measured at ambient pressure and room 
temperature (see Table [I]). This comparison allows us to 
analyze the effects of chemical pressure (i.e. anion substi- 
tution) on the electronic properties. For (TMTTF) 2 Br 
there are three additional Br bands crossing the lower 
organic band. 

We observe that the TMTTF bands vary only mod- 
estly with anion at a given temperature, particularly at 
the r and Z points. The largest difference in the band 
structure is seen along the X-V path, (X = (0.5,0,0), 
V = (0.5,0.5,0) in units of the reciprocal lattice vec- 
tors) where the indirect influence of the anion is most 
prominent. It is clear that in (TMTTF) 2 Br strong mix- 
ing with the Br bands distorts the TMTTF bands around 
the avoided crossings. Away from the avoided crossings 
the bands are similar to the TMTTF bands observed for 
the other salts. It is worth noting that (TMTTF) 2 Br is 
the only salt studied here with easily accessible metallic 
and superconducting statesP^ 

In Fig. [6] we show the real-space network of hop- 
ping terms (see Eq. [3]) between TMTTF molecules 
computed from the Wannier orbitals overlaps for 
(TMTTF) 2 AsF 6 (Fig.|6j(a)) and (TMTTF) 2 C10 4 (Fig. [6] 
(b)). The strength of the hopping is linearly encoded 
into the bond diameter. This figure shows that these 
materials have a preferred hopping direction (the di- 
rection with the thickest bonds), forming one dimen- 
sional chains. The ratio of inter-chain to intra-chain 
hopping strengths is smaller in (TMTTF) 2 AsFg than in 
(TMTTF) 2 C10 4 which indicates that (TMTTF) 2 AsF 6 is 
more one-dimensional. 
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X=SbF 6 , 7=300 K — 
X = Br, 7=300 K 



X= PF 6 , 7=300K — 
X=AsF 6 , 7=300K — 
X = BF d , 7=300 K 



FIG. 5. Band structures in the energy window [-0.8eV,0.2eV] 
of the TMTTF salts with crystal structures measured at am- 
bient pressure and room temperature. In this energy range 
all of the materials shown here have two bands arising from 
TMTTF HOMO orbitals. The (TMTTF) 2 Br salt addition- 
ally has three anion bands within this window. The com- 
mon TMTTF bands differ in details; between F and X (corre- 
sponding to the in-chain direction) the bands are very similar. 
There is more variation in the inter-chain direction, indicating 
the differing degrees of inter-chain coupling. 



While all of the materials have strong intra-chain hop- 
ping terms (of order ~0. 15-0.25 eV)), the inter-chain hop- 
ping terms can vary by about an order of magnitude (see 
Table [TT| in the Appendix). The values of the intra-chain 
hopping parameters in our work are consistent with those 
found for similar systems in previous experimental and 
theoretical investigations!*™^! Missing in those previous 
studies is a thorough analysis of the intra-chain dimer- 
ization as well as the inter-chain hopping parameters. 

In Table [I] we quantify the electronic dimerization for 
the Fabre CT salts studied in this work analogously as 
we did for the structural dimerization, i.e. 



dplrr 2 



to — 1\ 



(4) 



where to (t\) is the hopping term corresponding to the 
smallest (second smallest) bond length. While we ob- 
serve a significant dependence on the nature of the an- 
ion, the structural and electronic dimerizations seem to 
be uncorrelated. This can be understood physically: 
the electronic dimerization is defined by hopping inte- 
grals whose magnitude depends on the orientation of the 
overlapping orbitals as well as on their separation. If 
the orientation is more favorable along the longer intra- 
chain bond, then the more distant overlap can be larger. 
This is the case for the structures of (TMTTF) 2 PF 6 
above P = 0.9 GPa, (TMTTF ) 2 BF 4 at T = 100 K, 
and (TMTTF) 2 Br at room temperature; a negative value 



of d e i ec in Table [T] indicates that the longer bond has a 
larger hopping strength. Focusing on the anions with oc- 
tahedral symmetry at room temperature, we observe that 
while the structural dimerization has a consistent trend 
downwards as the anion changes from (SbFg) - through 
(AsFg) - to (PFg)~ (chemical pressure, smaller volume) 
and then further downwards as pressure is applied, the 
electronic dimerization follows the opposite trend. We 
will discuss this behavior below. 



B. Temperature dependence of the structural and 
electronic properties 

We proceed now with the analysis of the tempera- 
ture dependence of the structural and electronic behavior 
of a few Fabre CT salts. This study is done by per- 
forming ground state DFT calculations for structures 
determined at different temperatures. The investiga- 
tion for (TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 is done 
by considering crystal structures obtained experimen- 
tally at T = 4 K and at T = 300 K. The investigation 
for (TMTTF) 2 SbF6 is done using crystal structures de- 
termined experimentally at temperatures between T = 
100 K and 300 K. 

In Fig. [7] we present the band structure of 
(TMTTF) 2 SbF 6 as a function of temperature. We ob- 
serve that as the temperature is decreased, the band 
width increases and the dispersion between X and V be- 
comes steeper; this indicates that the electronic structure 
becomes more two dimensional with decreasing temper- 
ature. This trend can be also observed in the behavior of 
the 2D tight-binding parameters (Fig. [8J), especially i 2 , 
ti and £5. 

In order to quantify the electronic dimensionality, we 
introduce a dimensionality parameter, D, by taking the 
ratio of the inter-chain hopping terms (t„) and intra- 
chain hopping terms (t|), 



D 



EJtj\ 



(5) 



We emphasize that this parameter is an estimate of a 
model dimensionality. The correlation between tempera- 
ture, dimensionality, and bandwidth is seen more clearly 
by using this parameter, as illustrated in Fig. [9] (a), where 
D for (TMTTF) 2 SbF6 increases with decreasing temper- 
ature. 

In Fig. [TO] we present the band structure for 
(TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 for the crystal 
structures at T = 4 K and at T = 300 K. The inter- 
chain (X-V path) dispersion increases with decreasing 
temperature. (TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 un- 
dergo spin-Peicrls transitions (at T = 11.4 K and 16.4 K, 
respectively),^ however there is no energy splitting at 
T = 4 K since the crystal structure is not tetramerized. 
Interestingly, the electronic and structural dimerizations 
in these systems (see Table [TJ are larger for the room 
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(a) (TMTTF) 2 AsF 6 (b) (TMTTF) 2 CI0 4 




FIG. 6. Visualization of the strength of the hopping between the sites of the tight-binding model, the centers of mass of the 
TMTTF molecules (gray spheres); shown for (a) (TMTTF) 2 AsF 6 and (b) (TMTTF) 2 C10 4 (both at room temperature). The 
diameter of the bonds is proportional to the tight binding parameter strength \t a \. \t a \ above 0.1 eV are each a different shade 
of blue, \t a \ between 0.1 eV and 0.01 eV are a shade of red/orange, while \t a \ less than 0.01 eV are a shade of green. See 
Table |TT| in the Appendix for t a values. 
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FIG. 7. Band structure of (TMTTF) 2 SbF 6 calculated from 
the crystal structures obtained at several temperatures be- 
tween T — 100 K and 300 K. At T = 140 K and 300 K, struc- 
tures from two different samples were used; the additional 
bands at those temperatures are plotted with dashed lines. 
As the temperature is decreased, the bandwidth increases and 
the dispersion between the X and V points becomes steeper; 
this indicates that the electronic structure becomes more two 
dimensional with decreasing temperature. 



temperature structures than for the structures measured 
at T = 4 K. 

We also investigated (TMTTF) 2 BF 4 . For this system, 
the electronic dimcrization changes sign since the elec- 
tronic dimers are on the more closely spaced TMTTFs 
(in terms of center of mass separation) in the room tem- 
perature structure, and on the more distant pair for the 
structure at 100 K. 




7 d(ty (A) 



FIG. 8. Evolution of tight binding parameters of 
(TMTTF) 2 SbFe with temperature. As temperature is low- 
ered to T — 100 K, the dominant hoppings to an d t 1 become 
nearly equal, making the TMTTF chain nearly isotropic. The 
sizable 2D couplings t 2 , t<t and t$ show a complicated tem- 
perature dependence, with £ 2 and ts changing sign and tn 
increasing considerably as temperature is lowered. 



C. Pressure dependence of structural and 
electronic properties 

Here we investigate a series of new experimental crystal 
structures of (TMTTF^PFg determined at room tern- 
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FIG. 9. Dimensionality versus bandwidth of (TMTTF) 2 SbF 6 
with temperature (diamonds) and (TMTTF)2PF6 under pres- 
sure (circles); the bandwidth is the energy difference between 
the highest and lowest energies in the TMTTF bands around 
the Fermi energy, and the dimensionality is defined by the 
ratio of the hopping integrals in the intra- and inter-chain 
directions (see Eq. |5|). We see the expected positive cor- 
relation between bandwidth and pressure (indicated by the 
arrow): as the pressure increases, so does the intermolec- 
ular hopping and therefore the bandwidth. We also see a 
strong positive correlation between pressure and dimension- 
ality. There is a negative correlation between temperature 
and dimensionality and bandwidth; increasing the tempera- 
ture has a similar effect to decreasing the pressure. The 140 
K and 300 K points for (TMTTF) 2 SbF 6 are averaged over 
the two structures available at those temperatures. 



perature under various pressures. Figure [IT] shows how 
the band structure evolves as a function of pressure. As 
the pressure is increased, the bandwidth increases, and 
the system becomes more two-dimensional (i.e. the dis- 
persion is enhanced along the path X — V). This is also 



apparent in the tight-binding parameters (Figure 12); all 
the parameters grow with pressure (increasing the band- 
width), but not all by the same proportion, changing the 
degree of two-dimensionality. This trend to higher di- 
mensionality has also been observed experimentally - op- 
tical experiments on (TMTTF^PFg under pressure show 
that the metallic conductivity (Drude spectral weight) 
changes very anisotropically; it increases quickly with 
pressure in the perpendicular direction, while in the in- 
chain direction there is very little change.^ This trend 
with pressure has also been seen in other similar systems 
((TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 ), and identif ied as 
a cross-over from a quasi-lD system to a 2D metal.^^ 




r X V 
X=AsF 6 , 7= 300 K — 
X=AsF 6 , 7=4K 



y r z 

X= PF 6 , T= 300 K — 
"6' 



X=PF R 



r=4K 



FIG. 10. Band structures of (TMTTF) 2 AsF 6 (red) at 
T — 4 K (dashed) and room temperature (solid line), 
(TMTTF) 2 PF 6 (blue) at T = 4 K (dashed) and room temper- 
ature (solid line). The black lines are the bands resulting from 
the model Hamiltonian Eq. (J2|, parametrized by the Wannier 
orbital overlaps for the 8 shortest hops. In the model used 
there is no dispersion in the k z direction since the 8 shortest 
hops are all in the same plane, i.e. we have a two dimensional 
model. It is clear from the DFT bands that the interplanar 
coupling is small, which is why the 2D model fits so well. 
Note that T = 4 K is below both the charge ordering and spin 
Peierls transitions of (TMTTF) 2 AsF 6 and (TMTTF) 2 PF 6 . 




FIG. 11. Band structure of (TMTTF) 2 PF 6 at various pres- 
sures. As the pressure is increased, the bandwidth increases, 
and the dispersion becomes steeper between the X and V 
points; the system becomes more two dimensional. These 
trends are made obvious in Figure [9] (b). 



Table[T]shows that between P = 0.9 GPa and 1.5 GPa, 
the TMTTF molecules become almost equally spaced 
(in terms of the centers of mass) since 6 s t rU ct ~ 0. At 
P=1.5 GPa, the larger inter-chain hopping is no longer 
associated to t Q , but to ti; to avoid a discontinuity, we 
make an exception to the numbering of t a with ascending 




3 5 7 o^(A) 

FIG. 12. Tight binding parameters of (TMTTF) 2 PF 6 at var- 
ious pressures as a function of distance. Note that as the 
pressure is increased, the trend is for the t's to become larger. 
In Figure [9] (b) we see that the increases are such that the 
system becomes more two-dimensional. 



distance and refer, for pressures 1.5 to 2.7 GPa, to the 
largest hopping as to even though it belongs to the second 
nearest neighbour distance. The inter-chain hoppings to 
and t\ do not become equal around 1.2 GPa because even 
when the centers of mass are equally spaced, the stag- 
gering of the molecules in the chain means that the two 
hopping integrals are not equivalent. At pressures above 
P = 1.3 GPa, (TMTTF) 2 PF 6 is known experimentally 
to become metallic, and at lo w tem peratures undergoes 
a spin-density wave transition! 31 * 36 ! 

Figure M (b) shows how the dimensionality and band- 
width of (TMTTF) 2 PF 6 varies with pressure. We ob- 
serve the expected trend of increasing bandwidth under 
pressure (forcing the TMTTF molecules closer together, 
increasing their intermolecular interactions). We also see 
that physical pressure changes the bandwidth more, for 
a given change in dimensionality. 



VI. MODEL CALCULATIONS 



A. Exact Diagonalization of an Extended Hubbard 
Model 



In the previous section we obtained the network of in- 
teractions relevant for the Fabre CT salts by means of 
DFT calculations. We proceed now with model calcula- 
tions in order to analyze the effect of correlations in these 
materials. 

Since some of the phases realized in these materials 




FIG. 13. Schematic representation of the TMTTF molecules 
in the conducting plane for the Fabre salts here studied, (a) 
Balls correspond to the TMTTF molecules and ellipses indi- 
cate the dimers. The site and dimer labeling is shown, (b) A 
sketch of a possible charge- and magnetically-ordered state; 
the size of the circles represent the hole density, and the ar- 
rows the net magnetism. 



are charge and spin ordered phases, we shall investigate 
charge and spin structure factors using a quarter (hole) 
filled extended Hubbard model, 

H=- M C L<> + C ] a C i°) 

+U ^2 n it n H + ^2 VijUiUj, (6) 

where the sum over is over the 8 shortest dis- 

tances between sites, Uj (see Eq. [3]) are the correspond- 
ing hopping integrals, cJ CT (c^) is the creation (annihila- 



U and Vij are the 



tion) operator of a hole on the i th site with spin a, and 
7ij = nif + n s | with n ia = c\ a 
on-site and the inter-site Coulomb interactions 

In this Section we present numerical results for the 
Hamiltonian Eq. (|6]) , using the electron hopping parame- 
ters (t a with a = 0, 7) as defined in Eq.[3]and obtained 
in Sec. V. We choose the on-site Coulomb interaction to 
be typical for this class of materials U = 4io ~ 1 eV 
(U is of the order of the bandwidth W ~ 1 eV).^ Since 
the arrangement of the molecules changes only slightly 
with pressure and temperature, we assume that the pri- 
mary changes to the intersite interaction are based on 
the distance between the sites. This allows us to reduce 
the number of free parameters in our model; we scale 



V a as a function of the distance r a , V a = Vo^-, where 
the index a corresponds to the label of the hopping pa- 
rameter between that pair of sites. The intersite inter- 
actions Vj and V3 are set to zero since we expect these 
terms to be strongly screened by the intermediate sites. 
Thus we only have two 'free' parameters remaining, U 
and Vq. The Coulomb interaction along the chain Vq has 
been estimated in a previous work as between 0.2U and 
0.6t/.^ Here we consider two cases, Vb = 0.5to (weak 
intersite Coulomb repulsion) and Vq = 2to (strong inter- 
site Coulomb repulsion), both using U = 4i - With this 
set of parameters, the ground state for a system of size 
N = 16 (4x4) sites with periodic boundary conditions 
(see Figure [l3| is fou nd usi ng exact diagonalization, as 
implemented in ALPS EMUD While similar methods have 
been applied to some members of this family of materi- 
als before, we note that calculating our parameters from 
Wannier orbitals allows us to have a more c omple te, re- 
alistic description of the inter-chain coupling! 32 ! 40 ! 

We compute dimer structure factors for charge and 
spin, 



1, j 

with nj = (n. 



iq-(r. ;-r/) 



H+l 



)/2, 



Afb(q) =^E< m ^) e<, ' (r,_ri) 
d 1, j 

with mi — (mi + mj+i)/2, 



(7) 



(8) 



1 

ure 



where /, J are the dimer indices with i = 2(7 — 1) (and 
and i + 1 are the site (monomer) indices, see Fig- 
Nd is the total number of dimers, 17 denotes 
the position of the I th dimer, nj is the charge differ- 
ence between the sites in the dimer, and m/ is the total 
magnetization of di mers with m, = — the local 
magnetizationP2E2l Note that Cjj quantifies the correla- 
tion between the charge polarization of dimers, while M]j 
measures the correlation between spins on dimers. 

For the U and V values considered here, Cu(q) has 
a maximum at q = (0, 0) that corresponds to a charge 
order as shown in Fig. 13 (b). Mfl(q) has a maximum at 
q = (7r, 7r) corresponding to a dimer antifcrromagnetic or- 
der in both the in-cha in a nd i nter-c hain directions (shown 
(b))B2HH 
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schematically in Fig. 

Here we investigate charge and spin structure factors 
for various Fabre CT structures at different temperatures 
and pressures. To minimize the effects of experimental 
variability, we focus our analysis on the sets of struc- 
tures synthesized and measured following the same pro- 
cedure; the series of (TMTTF) 2 PF 6 under pressure, and 
(TMTTF) 2 SbF 6 for various temperatures. 

Both (0, 0) charge order and (tt, it) spin order are 
slightly suppressed with decreasing temperature and in- 
creasing pressure (Figs. 14 and[l5|). This is in contrast to 
a previous work on a simpler model which showed differ- 
ent trends for the charge and magnetic orders!^ Charge 
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FIG. 14. Structure factors as a function of temperature for 
(TMTTF)2SbFs. As the temperature is decreased the dimer 
charge and magnetic orders are somewhat suppressed. Simul- 
taneously, the bandwidth and dimensionality increase. Thus 
decreasing temperature has the same effects as increasing 
pressure (see Figure 15 I; the charge order is strongly acti- 



vated by Vb, and the magnetic order is only weakly enhanced. 
The squares correspond to results with Vb = 2to while the 
triangles correspond to Vb = 0.5to- 



order is strongly activated by increasing the strength of 
the inter-site Coulomb interaction, Vq. The antiferro- 
magnetic correlation is relatively weakly enhanced by in- 
creasing Vq. The changes in the correlation functions 



shown in Figs. 14 and 15 seem to be strongly connected 
to the degree of electronic dimerization d e i ec t r0 nic which 
we show in Fig. |16| Lowering the temperature down to 
T = 100 K in (TMTTF) 2 SbF 6 continuously decreases 
the electronic dimerization in these structures and thus 
suppresses intra-dimer charge order as well as inter-dimer 
antiferromagnetic order by making the one-dimensional 
chains more isotropic. The same observation holds for 
the increase of pressure on (TMTTF^PFe structures. 

The charge- and magnetically-ordered states found for 
the Fabre CT salts with this model are consistent with 
the phase diagram for these materials. But we find no 
evidence of a phase transition as a function of pressure or 
intersite Vs in this model. Note that this model cannot 
capture a spin-Peierls transition since no magnetoelastic 
coupling has been considered in the Hamiltonian. 

In this analysis we have concentrated on the (0, 0) 
charge order and (n, n) spin order. However, in prin- 
ciple, there may be many kinds of charge order in these 
systems. For example, a maximum at (0, tt) would indi- 
cate a charge order that alternates in the b direction as 
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FIG. 15. Structure factors as a function of pressure for 
(TMTTF^PFe. As the pressure is increased the dimer charge 
and magnetic orders are somewhat suppressed. Simultane- 
ously, the bandwidth and dimensionality increase (seen in Fig- 
ure[9](b)). The charge order is strongly activated by the value 
of Vo, while the magnetic order is only weakly enhanced, and 
is even finite for Vo = (not shown). The circles correspond 
to results with Vo — 2to while the diamonds correspond to 
Vo = 0.5t . 



well as in the a direction. Within the realistic parameter 
range explored here, we only observe the type illustrated 
in Fig. [131(b). 



VII. DISCUSSION 

Summarizing our results, we observe that the domi- 
nant band structure parameters obtainedin our work are 
generally consistent with those published for similar ma- 
terials, however Ref. 03] finds quite different values for 
the electronic dimerization. Those authors compute the 
t values by constructing TMTTF HOMOs from an ex- 
tended Hiickel model (a tight-binding model for both the 
a- and 7r-bonding systems of a molecule) and calculate 
the overlaps between them. This method does not al- 
low for the charge reorganization and other effects in the 
crystal, which are better included by using the Wannier 
orbitals from DFT. 

Table U shows that while structural dimerization tends 
to decrease with increasing pressure (both chemical and 
physical), the electronic dimerization only shows such 
a trend with physical pressure; there is no clear trend 
in electronic dimerization versus chemical pressure. Un- 
der chemical pressure, many aspects of the molecular ar- 




P(GPa) 



FIG. 16. Electronic dimerization for (a) (TMTTF) 2 SbF 6 
as function of temperature and (b) for (TMTTF)2PF6 under 
pressure. 



rangement can change (such as spacing and staggering), 
and there is no guarantee that they change smoothly 
with any one parameter of the anion, such as vol- 
ume. This is clearest in the series of anions (SbFg) - , 
(AsFg)~, (PFe)~ at room temperature. As the anion 
size decreases (equivalent to increasing chemical pres- 
sure) the electronic dimerization shows no trend, while 
the structural dimerization decreases. It is important to 
note that the dimensionality (also computed from the 
electronic hopping parameters) shows a clear trend of 
increasing with increasing chemical and physical pres- 
sure; a trend that has been observed expe rimentally in 
(TMTTF) 2 PF 6 and (TMTTF) 2 AsF 6 pEH As the tem- 
perature of (TMTTF)2SbF6 is increased, it shows a clear 
increase in both structural and electronic dimerization. 

Our model calculations show that while charge order is 
strongly activated by the inter-site Coulomb interaction, 
V, the magnetic order is weakly enhanced. We also see a 
weak suppression of both kinds of order as the pressure 
is increased, and as the temperature is decreased. 

To reproduce and understand the full phase diagram of 
these strongly correlated materials, one needs estimates 
of the Coulomb parameters. It is well known that molec- 
ular Coulomb parameters are overestimates for organic 
crystals; within the crystal, the interactions are strongly 
screened BUSH There are several promising approaches to 
calculating the screened Co ulomb parameters, each with 
their own costs and benefits! 46 * 47 * This will be addressed 
in a future work. 
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VIII. CONCLUSIONS 

We have examined the structural and electronic prop- 
erties of a set of Fabre charge transfer salts with crystal 
structures measured at different temperatures and pres- 
sures. By considering ab initio density functional the- 
ory calculations we obtain a comparable set of physically 
meaningful electron hopping parameters. In these results 
we identify some general trends: the structural dimeriza- 
tion is higher for the room temperature systems, the elec- 
tronic dimerization decreases with increasing pressure, 
the systems are more two dimensional at lower temper- 
atures and higher pressures, and this change in dimen- 
sionality is reflected in the degree of order in our model 
Hamiltonian. With this set of parameters, one can sys- 
tematically investigate the differences between these ma- 
terials in a model Hamiltonian. 



It is possible that the variations seen in the electronic 
structure (such as the change in dimensionality) are re- 
sponsible for tuning the ground states through the vari- 
ous phases accessible in these materials. However, ideally 
one would like a similarly systematic set of many-body 
interaction parameters, as well as the one-body parame- 
ters given here. 

Appendix A: Tight binding models 

In Table |H] we list all tight binding parameters we ob- 
tained from Wannier function overlaps. 

REFERENCES 



1 D. Jerome, Science 252, 1509 (1991). 

2 H. Mori, J. Phys. Soc. Jpn. 75, 051003 (2006). 

3 S. Yasuzuka and K. Murata, Sci. Technol. Adv. Mater. 10, 
024307 (2009). 

4 T. Kawakami, T. Taniguchi, S. Nakano, Y. Kitagawa, and 
K. Yamaguchi, Polyhedron 22, 2051 (2003). 

5 W. Yu, F. Zhang, F. Zamborszky, B. Alavi, A. Baur, C. A. 
Merlic, and S. E. Brown, Phys. Rev. B 70, 121101 (2004). 

6 N. Doiron-Leyraud, P. Auban-Senzier, S. R. de Cotret, 
K. Bechgaard, D. Jerome, and L. Taillefer, Physica B 405, 
S265 (2010). 

7 K. C. Kandpal, I. Opahle, Y.-Z. Zhang, H. O. Jeschke, and 
R. Valentf, Phys. Rev. Lett. 103, 067004 (2009). 

8 H. O. Jeschke, M. de Souza, R. Valentf, R. S. Manna, 
M. Lang,and J. A. Schlueter, Phys. Rev. B 85, 035125 
(2012). 

9 Single crystals of (TMTTF) 2 PF 6 and (TMTTF) 2 SbF 6 
were grown electrochemically in single or double H-type 
glass cells at room temperature. A constant voltage of 
0.9-1.1 V was applied between platinum electrodes of ap- 
proximately 5x10 mm 2 and 15x20 mm 2 respectively, re- 
sulting in a starting current through the solution of about 
14 /J.A. Needle-shaped single crystals of several millimeters 
in length and less than a millimeter in width are ready to 
harvest in about 1-4 weeks for (TMTTF) 2 SbF 6 , and 5-10 
weeks for (TMTTF) 2 PF 6 . 

Structural investigations on (TMTTF) 2 A with X =PF 6 
and SbF6 under ambient conditions and down to 100 K 
were performed at Universitat Stuttgart using a Kappa 
CCD Bruker AXS diffractometer. Diffraction angles 9 of 
0.41°-28.28° were considered while irradiating with a wave- 
length of A=0. 71073 A. The obtained R-values vary from 
0.036 to 0.062. 

The pressure dependent x-ray diffraction data for 
(TMTTF) 2 PF 6 were collected at the ID09A beamline of 
the European Synchrotron Radiation Facility in Greno- 
ble using a diamond anvil cell. A wavelength of 0.413 A 
was irradiated with a diffraction angle 9 of about 25°. X- 
ray diffraction patterns were collected on an imaging plate 
MAR345 detector by rotating the crystal from —30° to 



+30° with 2° steps, and analyzed using the XDS packaged 
The refinement of the atomic positions was performed by 
the SHELX software! 21111 

Gaseous helium was used as the pressure transmitting 
medium; the pressure in the cell was determined in situ 
by the ruby luminescence method.^ Typical R-values are 
in the range of 0.07 - 0.08. 

10 W. Kabsch, J. Appl. Cryst. 26, 795 (1993). 

11 A. Pashkin, M. Dressel, S. G. Ebbinghaus, M. Hanfland, 
and C. A. Kuntscher, Synth. Met. 159, 2097 (2009). 

12 E. Rose, C. Loose, J. Kortus, A. Pashkin, C. A. Kuntscher, 
S. G. Ebbinghaus, M. Hanfland, F. Lissner, T. Schleid, and 
M. Dressel, J. Phys.: Condens. Matter 25, 014006 (2013). 

13 H. K. Mao, J. Xu, and P. M. Bell, J. Geophys. Res. 91, 
4673 (1986). 

14 K. Koepernik and H. Eschrig, Phys. Rev. B 59, 1743 
(1999). 

15 J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

16 G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 

17 G. Kresse and J. Furthmiiller, Comput. Mat. Sci. 6, 15 
(1996). 

18 P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

19 G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

20 S. Grimme, J. Comput. Chem. 27, 1787 (2006). 

21 T. Granier, S. Gallois, L. Ducasse, A. Fritsch, and A. Fil- 
hol, Synth. Met. 24, 343 (1988). 

22 B. Liautard, S. Peytavin, G. Brun, and M. Maurin, Cryst. 
Struct. Commun. 11, 1841 (1982). 

23 J. L. Galigne, B. Liautard, S. Peytavin, G. Brun, J. M. 
Fabre, E. Torreilles, and L. Giral, Acta. Cryst. B 34, 620 
(1978). 

24 B. Liautard, S. Peytavin, G. Brun, D. Chasseau, J. M. 
Fabre, and L. Giral, Acta. Cryst. C 40, 1023 (1984). 

25 J. L. Galigne, B. Liautard, S. Peytavin, G. Brun, M. Mau- 
rin, J. M. Fabre, E. Torreilles, and L. Giral, Acta. Cryst. 
B 35, 1129 (1979). 

26 D. Pedron, R. Bozio, M. Meneghetti, and C. Pecile, Phys. 
Rev. B 49, 10893 (1994). 

27 P. M. Grant, Phys. Rev. B 26, 6888 (1982). 



12 



Anion 


A* 


to ti 


t 2 


tz £4 


ts 


t 6 


tr 


optimized 




SbF 6 100 K 


-0.2516 


0.1823 0.1747 


-0.0344 


-0.0030 0.0310 


-0.0104 


0.0018 


-0.0027 


yes 




SbF 6 140 K (1) 


-0.2459 


0.1810 0.1692 


-0.0285 


-0.0029 0.0314 


-0.0061 


0.0009 


-0.0027 


yes 




SbF 6 140 K (2) 


-0.2486 


0.1853 0.1687 


-0.0279 


-0.0029 0.0318 


-0.0059 


0.0008 


-0.0030 


yes 




SbF 6 180 K 


-0.2449 


0.1853 0.1652 


-0.0193 


-0.0029 0.0283 


0.0010 


0.0001 


-0.0020 


yes 




SbF 6 200 K 


-0.2417 


0.1857 0.1612 


-0.0138 


-0.0029 0.0259 


0.0053 


-0.0004 


-0.0013 


yes 




SbF 6 300 K (1) 


-0.2379 


0.1947 0.1470 


0.0023 


-0.0032 0.0176 


0.0146 


< ±10~ 4 


-0.0010 


yes 




SbF 6 300 K (2) 


-0.2342 


0.1925 0.1426 


0.0043 


-0.0034 0.0140 


0.0129 


0.0012 


-0.0007 


yes 




AsF 6 4K 


-0.2617 


0.1943 0.1759 


-0.0380 


-0.0038 0.0358 


-0.0141 


0.0015 


-0.0051 


no 




PF 6 4K 


-0.2539 


0.1912 0.1686 


-0.0333 


-0.0038 0.0367 


-0.0104 


0.0008 


-0.0052 


no 




AsF 6 300 K 


-0.2292 


0.1751 0.1568 


-0.0076 


-0.0047 0.0188 


0.0244 


0.0010 


-0.0021 


yes 




PF 6 300 K 


-0.2477 


0.1976 0.1569 


-0.0025 


-0.0141 0.0030 


0.0312 


-0.0003 


-0.0035 


yes 


PF e 


300 K, 0.3 GPa 


-0.2203 


0.1981 0.1093 


0.0059 


-0.0056 0.0278 


-0.0094 


0.0015 


-0.0063 


yes 


PF e 


300 K, 0.6 GPa 


-0.2280 


0.2065 0.1118 


0.0108 


-0.0062 0.0291 


-0.0090 


0.0009 


-0.0066 


yes 


PF e 


300 K, 0.9 GPa 


-0.2513 


0.2193 0.1348 


0.0019 


-0.0065 0.0349 


-0.0085 


0.0003 


-0.0085 


yes 


PF e 


300 K, 1.5 GPa 


-0.2550 


0.2207 0.1390 


0.0015 


-0.0080 0.0378 


-0.0136 


0.0001 


-0.0086 


yes 


PF e 


300 K, 2.0 GPa 


-0.2757 


0.2333 0.1561 


-0.0041 


-0.0092 0.0408 


-0.0159 < ±10~ 4 


-0.0093 


yes 


PF e 


300 K, 2.7 GPa 


-0.3081 


0.2398 0.1996 


-0.0199 


-0.0092 0.0513 


-0.0154 


-0.0015 


-0.0099 


yes 




Br 300 K 


-0.2215 


0.1719 0.1422 


-0.0270 


-0.0048 0.0311 


0.0025 


-0.0007 


-0.0005 


yes 




C10 4 300 K (A) 


-0.2195 


0.2017 0.1067 


-0.0022 


-0.0042 0.0313 


-0.0136 


0.0010 


-0.0081 


no 




CIO4 300 K (B) 


-0.2289 


0.2017 0.1067 


-0.0022 


-0.0050 0.0320 


-0.0136 


0.0011 


-0.0081 


no 




BF 4 100 K (A) 


-0.2435 


0.1735 0.1644 


-0.0151 


-0.0063 0.0296 


-0.0240 


0.0004 


0.0210 


no 




BF4 100 K (B) 


-0.2353 


0.1735 0.1644 


-0.0151 


-0.0062 0.0300 


-0.0240 


0.0015 


0.0210 


no 




BF4 300 K (A) 


-0.2618 


0.2057 0.1466 


-0.0209 


-0.0010 0.0381 


-0.0220 


-0.0013 


-0.0032 


no 




BF4 300 K (B) 


-0.2413 


0.2057 0.1466 


-0.0209 


-0.0007 0.0373 


-0.0220 


-0.0013 


-0.0032 


no 



TABLE II. tot values determined from the Wannier orbitals for all TMTTF structures investigated here (energies in eV). The 
t a are numbered from shortest to longest bond (defined by distance between the centres of mass of the TMTTF molecules) , 
except where noted below. It is clear that the intra-chain terms, to and t\, are the dominant hopping terms. Note that for 
(TMTTF)2PF6 above 0.9 GPa, for Br at room temperature, and for BF4 at 100 K the dominant in-chain t is the longer one; 
we have swapped the labels for these materials such that to remains the strongest in-chain coupling. The labels (A) and (B) 
refer to the two inequivalent TMTTFs in each unit cell (in the absence of inversion symmetry) . 
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